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Abstract 

Once formed in a supernova explosion, a neutron star cools rapidly via neutrino emission 
during the first 10 4 -10 5 yr of its life-time. Here we compute the axion emission rate from 
baryonic components of a star at temperatures below their respective critical temperatures 
T c for normal-superfluid phase transition. The axion production is driven by a charge neutral 
weak process, associated with Cooper pair breaking and recombination. The requirement 
that the axion cooling does not overshadow the neutrino cooling puts a lower bound on the 
axion decay constant f a > 6 x lO 9 !^ 1 GeV, with T c g = T c /10 9 K. This translates into a 
upper bound on the axion mass m a < 10~ 3 T c g eV. 
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1. Introduction 

CP-violation in the strong sector of the Standard Model arises due to a topological 
interaction term in the QCD Lagrangian 

where F^ u = <9 M v4„ — <9„v4 M + gf abc A^f,A uc is the gluon field strength tensor, g is the strong 
coupling constant, F* = e^ u \ p F xpa /2, f abc are the structure constants of SU(3) group, 9 
is the parameter which parametrizes the non-perturbative vacuum states of QCD \9) = 
5^ n exp(— in6)\n), where n is the winding number characterizing each distinct state, which 
is not connected to another by any gauge transformation [lj]. The QCD action changes by 
27r under the shift 9 — > 9 + 2tt, i.e., 9 is a periodic function with a period of 2tt. In presence 
of quarks the physical parameter is not 9, but 

9 = 9 + argdetm,, (2) 

where m q is the matrix of quark masses. Experimentally, the upper bound on the value 
of this parameter is 9 < 10~ 10 , which is based on the measurements of the electric dipole 
moment of neutron d n < 6.3 • 10~ 26 e cm 0]. The smallness of 9 is the strong CP problem: 

Preprint submitted to Nucl. Phys. A, vol. 897, Jan. 2013, pp. 62-69 31 May 2012 



the Standard Model does not provide any explanation on why this number should not be of 
order unity. 

An elegant solution to the strong CP-problem is provided by the Peccei-Quinn mecha- 
nism This solution amounts to introducing a global U(1)pq symmetry, which adds 
an additional anomaly term to the QCD action proportional to the axion field a. This term 
acts as a potential for the axion field and gives rise to an expectation value of the axion 
field (a) ~ —9. The physical axion field is then a — (a), so that the undesirable 9 term in 
the action is replaced by the physical axion field. The axion is the Nambu-Goldstone boson 
of the Peccei-Quinn U(1)pq symmetry breaking 0,0], and its effective Lagrangian has the 
form 

& a = ~\d^a + J2U(0 M a, V), (3) 

where the second term describes the coupling of the axion to fermion fields (ip) of the 
Standard Model. 

There are ongoing experimental searches for the axion and the cosmology and astro- 
physics provide strong complementary constraints. Because axions can be effectively pro- 
duced in the interiors of stars they act as an additional sink of energy. The requirements 
that the energy loss from a star is consistent with the astrophysical observations place lower 
bounds on the coupling of axions to the Standard Model particles, and hence on the Peccei- 
Quinn symmetry breaking scale j6[ 0]. The latter limit translates into an upper limit on 
the axion mass. Such arguments have been applied to the physics of supernova explosions 
[8 12] and white dwarfs 13]. In the case of supernova explosions the dominant energy loss 
process is the emission of an axion in the nucleon (n) bremsstrahlung n + n — > n + n + a. 
The same process was considered earlier by Iwamoto as a cooling mechanism for mature 



neutron stars, i.e., neutron stars with core temperature in the range 10 8 — 10 9 K [14J. The 



implications of the axion emission by the modified nucleon bremsstrahlung, as calculated 



in Ref. 14|, on the cooling of neutron stars were briefly discussed in Ref. 15]. However 



it is now well established through cooling simulations of compact stars 16Ml8l] . that their 
neutrino cooling era, which spans the time period t < 10 4 — 10 5 yr after birth, is strongly 
affected by the neutrino emission from its superfluid phases due to the process of neutrino 
emission by Cooper-pair breaking 19M24 . 

In this article we compute the rate at which the superfluid phases of a neutron star lose 
their energy by axion emission via the processes of Cooper pair breaking and recombination. 
This work concentrates on the baryonic interiors of compact stars and considers for the 
sake of simplicity S'-wave superfluids. The inner crusts and the baryonic core of a neutron 
star features iso-triplet spin-0 S- wave superfluids; in addition the core may contain spin-1, 
P-wave neutron superfluid 25M27I] and at high densities, iso-singlet, spin-1, neutron-proton 
.D-wave superfluid [28|. The extension of the present work to P-wave and P-wave superfluids 
is straightforward. We use natural units, h — c — fee = 1- 



2 



2. Currents, matrix elements and emissivity 

The coupling of axion to baryonic fields is described by the following interaction La- 
grangian 



&int = (4) 



1 

fa 

where f a is the axion decay constant, the baryon and axion currents are given by 



B» = Cj 7 ^, L, = d^, (5) 

where C a are model dependent coupling constants of order of unity. In the case of a mul- 
ticomponent baryonic system the baryon current contains a sum over all components. The 
squared matrix element for the process of axion emission is then given by 

\M a \ 2 = i/- 2 (S^)(L M Lt). (6) 

The energy radiated per unit time in axions (axion emissivity) is given by the phase-space 
integral over the probability of the process of emission 

f d 3 q 

e a = -f~ 2 J ug(u)q^q u Imn^(q), (7) 

where q and u are the axion momentum and energy. Here we defined the polarization tensor 
of baryonic matter 

n i 

where the i sum is over the four-momenta of the baryons. Upon carrying out the angular 
integral in Eq. (J7J) we write the emissivity in terms of a one- dimensional integral 



fa 



47T 2 



d\q\ q g(u)K a (q), (9) 



where the contraction of the axion currents with the baryonic polarization tensor is given 
by 

« a (g) = g^JmlTf (g). (10) 

So far the expression for the axion emissivity is completely general; we will need to compute 
the polarization tensor of baryonic matter for the process of interest. 



3. Polarization tensor of superfluid baryon matter 

At sufficiently low densities and temperatures baryonic matter forms a 1 Sq pair con- 
densate. In compact stars this is the case for all baryons except neutrons, which may form 
P-wave superfluid at densities at and above the saturation density. To describe the response 
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Figure 1: The two diagrams contributing to the polarization tensor of baryonic matter, which defines 
the axion emissivity. The "normal" baryon propagators for particles (holes) are shown by single-arrowed 
lines directed from left to right (right to left). The double arrowed lines correspond to the "anomalous" 
propagators F (two incoming arrows) and F + (two outgoing arrows) . The horizontal dashed lines represent 
the axion a. 



of baryonic matter to the axion field we use the methods developed for the description of neu 



trino interactions in Refs. |22l. l24j (see also |23l. l29| ) . The imaginary-time momentum-space 



correlators are given by the 2x2 Nambu-Gor'kov matrix 

y <7j<7 >{iu n ,p) = ~ * . (11) 

V F av'( lUJ n,P) Gj a ,(lU n ,p) J 

The elements of the matrix are time-ordered correlators of the baryon field ips and i/j b ; in 
the frequency-momentum domain these are given by 

/ u 2 v 2 \ 

G aa '(iu n ,p) = 8 aa ,(, g— + . g , (12) 

F aa >(iu n ,p) = -ia y u p v p (- — — ], (13) 

where F^ o . l (iu n ,p) = F ac i (icu n , p) , uj n = {2n + l)nT is the fermionic Matsubara frequency, 
cr y is the y component of the Pauli-vector in spin space, u 2 = (1/2) (1 + ^p/e p ) and u 2 + 
v 2 = 1 Bp = + /S. 2 is the single particle energy in the paired state, with A p being 
the (generally momentum- and energy-dependent) gap in the quasiparticle spectrum and 
£p = vp{j> — pf) is the single-particle spectrum in the normal state, where vp and pf 
are the (effective) Fermi velocity and momentum. The propagator for the holes is de- 
fined as G^ al (iu n ,p) = G aa r(—i(jj n ,—p). For an 5-wave condensate we have G aa /(iu n ,p) = 
5 aa iG(iuj n ,p), F aa i(ioj n ,p) = —ia y F(iuj n ,p), etc. The polarization tensor of a superfluid ob- 
tains contributions from four distinct diagrams that can be formed from the normal and 
anomalous propagators with four distinct effective vertices. However, for the axial vector 
perturbations the vertices are not renormalized in the medium and, therefore, one proceeds 
with the bare vertices, in which case the number of the distinct contributions to the polar- 
ization tensor reduces to a sum of two admissible bare loops (see Fig. [T]) 

A T (q) = U GG (q)-niFF(q), (14) 
where T = ±1 is the time reversal operator and 



n GG (g) = T I t^J2 g ^ g (p + ^ = ^t [dx(G*G) : 



(15) 



where T is the temperature, p = {ip n ,P) with p n being the fermionic Matsubara frequency, 
i/(0) = m*pF/^ 2 is the density of states of neutrons with m* being their effective mass, 
which may include the wave function renormalization (the so-called E-mass), x is the cosine 
of the angle formed by the vectors q and p. The convolution is defined as 

/oo 
d^G(p)G(p + q), (16) 
-oo ■ 

with a similar expression for (F * F). Carrying out the Matsubara sums appearing in the 
convolutions we find 



A ± (q) = -2r/(0) 



x5 (1±1 
+ „ 



1-xS 



(F*F+), (17) 



where (...) = (1/2) J\ dx(...), x = q ■ p and 5 = \q\vp/oj is a small parameter of the 
theory. To leading order in 5 parameter the imaginary part of the (F * F) convolution can 
be computed analytically 



, ttA 2 6(u-2A) , / u \ 

whereas the real part follows from the dispersion relation 

Re(F*F + ) = - / lm{F*F + ) {u'). (19) 

n J-oo u-u' 

For the contraction fllOl) we obtain 

K a (q) = uHmU \q)-q i ulmU l0 (q)-uq j lmU 0l ( y q) + ^-lmn ll ( y q), (20) 

o 

where, to the accuracy we are working, lmU°^(q) = A + v%, Imn"(g) = 3A~(q), qjlxall^ (q) = 
A~ qvp and gjmil l0 (g) = A + qvF- Here A ± = A ± (q-p) are the first moments of A ± integrals 
with respect to the cosine of the angle formed by the vectors q and p. The leading order in 
S expansion of the contraction gives 



K a (q) = -2i/(0) 



0O 2 l'l ( 1 + " 



5 2 \ 2u 2 5 2 q"5 2 



Im(F * F H 



Aui0) fv F lm(F*F + ) + O(6 4 ). (21) 

d 

Eq. fT2TT) completes our evaluation of the baryonic polarization tensor and its contraction 
with axionic current, which is the key input for the computation of the axion emissivity. 
Note that to 0(5 4 ) accuracy the (F * F + ) convolution appears at its lowest (5 = 0) order. 
More accurate evaluation is not required, because of the major uncertainty in the coupling 
strength of axions to the Standard Model fermions as well as other large uncertainties in the 
physics of cooling neutron stars. 
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4. Axion emissivity 



The axion emissivity is obtained on substituting Eqs. f|T8|) and (l2Tj) in Eq. (Q 

e a = f f- 2 u(0)vlT"l aj (22) 



where 

/ a = z 5 I dy 1— te (zyf , 



„3 



i Vv 2 - 1 

z = A(T)/T and = [1 + exp(x)] -1 is the Fermi distribution function. The T 5 scaling 

of the emissivity is understood as follows. The integration over the phase space of neutrons 
carries a power of T, since for degenerate neutrons the phase-space integrals are confined 
to a narrow strip around the Fermi surface of thickness T. The axion is emitted thermally 
and being relativistic contributes a factor T 3 to the emissivity. The one power of T from 
the energy of the axion and the inverse one power of T from the energy conserving delta 
function cancel. The transition matrix element is proportional to the combinations of u p 
and v p amplitudes, which are dimensionless, but contain implicit temperature dependence 
due to the temperature dependence of the gap function. This dependence is not manifest in 
Eq. fl22|) . i.e., was absorbed in the definition of the integral I a . Thus, the explicit temperature 
dependence of the axion emission rate Eq. fj22|) is T 5 . In the cgs units the axion emissivity 
Eq. (J22J) is 

21 /10 10 GeV\ 2 fm*\ 2 sv F \Z ( T \ 5 3 , 

e a = 1.06x10-^— j (-] (f) (—) /.erg cm- b- (24) 

where two powers ofvp/c arise from the small momentum transfer expansion and one power 
- from the density of states. At temperatures of order the critical temperature T c ~ 10 9 K the 
superfluid cools primarily by emission of neutrinos via the pair-breaking processes driven by 
the axial-vector currents (we continue to assume that potential fast cooling via direct Urea 
processes is prohibited). The emissivity of this processes in the case of 1 5 , -wave superfluid 



is given by [19|, |20|, |23| 



( A v(0)v F 2 T 7 h, (25) 



15tt 3 

where Gp is the weak Fermi coupling constant, g A = 1-25 is the axial-vector current coupling 
constant, (a = 6/7 and 

I v = z 7 / dy JL— fr {zyf . (26) 

We now require that the axion luminosity does not exceed the neutrino luminosity, i.e., 

e a lOvr 2 I a 



e v flG 2 F g 2 A (,A h 



< 1. (27) 
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Substituting the the free-space value of the axial vector coupling g& = 1.25 and introducing 
r(z) = z 2 (I a /I u ) we transform Eq. ( |27|) 



59.2 



r(z). 



(2? 



Not far from the critical temperature A(T) ~ 3.06T c -^/l — T/T c , which translates into z = 
3.06 t _1 -\/l — i, where £ = T/T c . Numerical evaluations of the integrals provides the following 
values r(0.5) = 0.07, r(l) = 0.26, r(2) = 0.6 and asymptotically r{z) — > 1 for z ^> 1. 
Substituting the value of the Fermi coupling constant Gp = 1.166 x 10 -5 GeV -2 in Eq. 
and noting that r(z) < 1 , we finally obtain 



/ a > 5.92 x 10 9 GeV 



0.1 MeV 



A(T) j 

which translates into an upper bound on the axion mass 

10 10 GeV\ 



0.62 x 10^eV 



fa 



< 1.05 x 10" d eV 



A(T) 
0.1 MeV 



(29) 



(30) 



The bound (1291) can be written in terms of the critical temperature by noting that A(T) ~ T c 
in the temperature range 0.5 < t < 1 of most interest. 



5. Discussion 

The neutrino cooling era of compact stars, which spans the time-period t < 10 4 — 10 5 
yr after their birth in supernova explosions is a sensitive probe of the particle physics of 
their interiors. If one assumes that there are no rapid channels of cooling in neutron stars, 
i.e., deconfined quarks, above Urea threshold fractions of protons or hyperons (all of which 
lead to a rapid Urea cooling), then neutron stars cool primarily by neutrino emission in 
Cooper pair-breaking processes in baryonic superfluids. Here we have shown that if axions 
exist in Nature, the neutron stars must cool via axion emission in Cooper pair-breaking 
processes, whose axion emission rate scales as T 5 . This scaling differs from the T 7 scaling 
of the counterpart neutrino processes. The difference arises from the different phase spaces 
required for an axion and a pair of neutrinos and is independent of the baryonic polarization 
tensor. Note also that the rate of axion emission from a P-wave superfluid will differ from the 
S'-wave rate, derived above, by a factor 0(1) and, therefore, will not change quantitatively 
the obtained bounds on the axion parameters. 

If the Standard Model physics provides a consistent explanation of the data on the cooling 
of neutron stars, one can place a lower bound on the axion decay constant (the breaking 
scale of the Peccei-Quinn symmetry). Our calculations show that this bound is given by 
f a > 5.92 x 10 9 T~g GeV, where T c9 is the magnitude of the critical temperature in units of 
10 9 K. This translates into an upper bound on the axion mass m a < 10~ 3 T c9 eV. Similar 



bounds were obtained previously by Iwamoto 14J (//10 10 GeV > 0.3) from comparison of the 



rates of axion bremsstrahlung and modified Urea neutrino emission by mature neutron stars, 
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and by Umeda et al [15( (//10 10 GeV > 0.1 —0.2) from fits of cooling simulations to the PSR 
0656+14 dataQ Our lower bound on f a is somewhat larger than the one that follows form the 
requirement that the axions do not "drain" too much energy from supernova process so that 
it fails 0,0-12]. Furthermore, our bound sensitively depends on the pairing gap in baryonic 
superfluids, whose magnitude and density dependence are not well-known. The physical 
implications of the bound ( 12 9 p can be fully explored with numerical simulations of neutron 
star cooling. Targeted fits to a specific object exhibiting slow cooling would be more useful 
than fits to the entire population of neutron stars with measured X-ray fluxes. Examples, 



of such fits were carried out, e. g., in the case of the neutron star in CAS A [30|, |31j. The 
accuracy of the predictions will be limited by the uncertainties in the physical input in 
cooling simulations and uncertainties inherent to the interpretation of the data. For stars in 
the neutrino cooling era the dissipative heating processes are unimportant. The potential 



sources of uncertainty are well documented in the literature 16Ml8| and include (a) the rate 



of neutrino/axion emission, in particular, its dependence on the magnitude of the pairing 
gap; (b) the composition of the surface layers; (c) the influence of the star's 5-field on the 
photon and neutrino emission processes; (d) the estimate of the age of any given object. We 
anticipate that the error bars on the bounds should remain within a factor of few and the 
accuracy of the bounds can be improved with the help of numerical simulations of cooling 
compact stars. 
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